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Abstract 

In this paper, in the framework of nonrelativistic QCD we study the hght hadron (LH) decays 
of the spin-triplet (5=1) D-wave heavy quarkonia. The short-distance coefficients of all Fock 
states in the ^Dj{J = 1,2,3) quarkonia including the D-wave color singlet, P-wave color octet, 
and S'-wave color singlet and color octet are calculated perturbatively at order. The operator 
evolution equations of the four-fermion operators are also derived and are used to estimate the 
numerical values of the long-distance matrix elements. We find that for the cc system, the LH 
decay widths of ^p{l^Dj) predicted by nonrelativistic QCD is about 2 ~ 3 times larger than the 
phenomenological potential model results, while for the bb system the two theoretical estimations 
of r(T(l'^Dj) — )• LH) are in coincidence with each other. Our predictions for tpll^Dj) LH decay 
widths are r{^P{l^DJ) LH) = (435, 50, 172)keV for J=l,2,3; and for T{l^Dj), r{T{l^Dj) 
LH) = (6.91, 0.75, 2.75)keV for J=l,2,3. 
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I. INTRODUCTION 



The production, decay, and mass spectrum of heavy quarkonium have been interesting 
topics since the first charmonium state J/ip was discovered in 1974. Because of their large 
mass scales and nonrelativistic nature, heavy quarkonia are good probes to study and un- 
derstand quantum chromodynamics (QCD) from both perturbative and nonperturbative 
aspects. In fact, one of the earliest applications of QCD is to calculate the inclusive decay 
rates of heavy quarkonia. In early times, it was assumed that such a decay process can 
proceed through two steps. First, the heavy quarkonium transforms into a free QQ pair, 
which is a long-distance nonperturbative effect. Then, the heavy-quark pair annihilates into 
light hadrons (LH) through gluons, which can be calculated perturbatively. In the nonrel- 
ativistic limit, the long-distance part is related to the QQ Schrodinger wave functions or 
their derivatives at the origin. In this picture, the free QQ are in color singlet and have 
the same quantum numbers J^*" as the bound state heavy quarkonium. This is referred 
to as the "color-singlet model". Explicit calculations at next-to-leading order (NLO) in 
for S-wave quarkonium decays support the color^singlet model factorization formula. But 



it breaks down in the calculations of P-wave 



2| and D-wave 3|, |j| heavy quarkonium LH 



decays at order, where infrared divergences appear. Phenomenologically, these infrared 
divergences are regularized by the binding energy of QQ bound states. 

In Ref. |5j], Bodwin, Braaten, and Lepage first introduced the color-octet matrix elements 
to absorb the infrared logarithms, then they developed nonrelativistic QCD (NRQCD) ef- 
fective field theoryj^, based on which the inclusive decay rate of heavy quarkonium can 
be given by a rigorous factorization formula, and calculated in a systematic way by double 
expansion of a^, the coupling constant of QCD, and f , the typical velocity of heavy quarks 
in the heavy quarkonium. In their formula, heavy quarkonium is treated as a superposition 
state of \QQ), \QQg), \QQgg), and other higher order Fock states, rather than the \QQ) 
color-singlet state only. The contribution of each Fock state is organized in powers of f^, 
and can be written as a product of the long-distance matrix element and the corresponding 
short- distance coefficient. Huang and Chaofr] first got the infrared finite LH decay width 
of the spin singlet P-wave state, with QCD radiative corrections in the framework of 
NRQCD. The decay widths of Xcj LH were calculated to order in Refs.j^, 9]. Com- 
plete and detailed results of color-singlet and octet-short distance coefficients of S'-wave and 
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P-wave spin-triplet states were given in Ref . jlo| . In Refs.|lOl. the authors also explained 
why the infrared divergences disappear in the NRQCD factorization approach. 

NRQCD is now a widely accepted effective field theory for heavy quarkonium. In the 
framework of NRQCD, lots of theoretical work has been done to study 5*- and P-wave 
quarkonium decays, and some significant successes have been achieved (for a review see 
Ref. [l^ ) ■ Recently, the order results of S- and P-wave heavy quarkonium inclusive 



hadronic decays were obtained by Brambilla et al.[13| However, compared with S- and P- 
wave heavy quarkonium, little work has been done on D-^NdNe states in NRQCD. In this 
paper, we will calculate the LH decay widths of spin-triplet D-wave states ^Dj (J=l,2,3). 
Here the subprocesses ^S*!^'^^ LH and ^Pj^^ LH at leading-order (LO) in are all 
included, in which the short-distance coefficients are calculated in a different way and in 
agreement with the results in Ref. jl^. As the main part of these D-wave quarkonium decays, 
the infrared safe short-distance coefficients of ^d't' — t- LH are first obtained in this paper. 

n 

We use the covariant projection method, which was first introduced inp^ and generalized 



in Refs. jiol. fisjl. to do the perturbative calculations. At LO in f^, the long-distance matrix 
elements of color-singlet D-wave four-fermion operators are related to the wave function's 
second derivative at the origin. And the matrix elements of the S-wave and P-wave octet 
and of S'-wave singlet four-fermion operators could be studied in lattice simulations, or 
determined by fitting experimental data, or roughly estimated through velocity scaling rules. 
To give numerical predictions, here we use operator evolution equations to estimate the 
values of the matrix elements. The rest of our paper is organized as follows: In Sec. II, we 
briefly introduce the NRQCD effective field theory and give the general formulas used in 
D dimension. In Sec. HI, all the subprocesses will be calculated to order. And after 
matching the full QCD results with the NLO NRQCD ones, the infrared safe short- distance 
coefficients as well as the operator evolution equations are obtained in Sec. IV. In Sec.V, 
we will discuss the numerical results and their phenomenological applications to cc and bb 
systems. 

II. GENERAL FORMULAS 

The Lagrangian for NRQCD isj^: 

^NRQCD = flight + ^heavy + 5C, (l) 
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where Cught includes gauge field and light quark parts, and Cheavy is the nonrelativistic 
Lagrangian for the heavy quarks and antiquarks: 

= i-HiD. + + x'm - ^)x, (2) 

where ip ^ind ^I's the Pauli spinor fields that annihilate heavy quark and antiquark, 
respectively, and Dt and D are the time and space components of the gauge-covariant 
derivative D^. 6C describes the relativistic effects. The leading-order f ^ corrections are the 
bilinear terms: 

+ ^[^^D -gB-E- gB)^ + x^D -gE-B- gB)x] 

+ X gE-gEx iD) ■ ct^ + xH^D x gE - gE x iD) • ax] 

+ ^[^\9B ■ ct)^ - X^gB ■ ct)x], (3) 



where Cj could be obtained by using the Foldy-Wouthuysen-Tani transformation |16j . di- 
agonalizing the Dirac theory so as to decouple the heavy-quark and antiquark degrees of 
freedom. To reproduce the full QCD Lagrangian and to describe the annihilation of the 
heavy-quark pair, four-fermion local operator terms are also needed: 

SCi-fermion = ^H^-i ^nil^A), (4) 

n V 

where Ori(yUA) is the local four-fermion operator with the general form {Tp'^lCix){x^^'i'P)) 
and fni^^h) is the Wilson coefficient. In effective theory both operators and coefficients 
are dependent on the factorization scale /iA, but their combinations cancel the dependence. 
With the help of the optical theorem the inclusive decay rate of a heavy quarkonium state 
H could be expressed as: 

V{H ^ LH) = 2\m{H\5C^.f,rmion\H) = V ^^^^^|^(/7|a(/iA)|i^). (5) 

In principle, we need infinite terms to give theoretical predictions, but in practice only a finite 
number of these terms are needed to give an order of f ^ result, since the long-distance matrix 
elements can be ordered in powers of f ^ by applying the velocity scaling rules summarized 
in Ref. 6|. And the short-distance coefficients (Wilson coefficients), defined in the matching 



condition below, can be calculated perturbatively as a perturbation series in QCD coupling 
constant 

AiQQ^QQ) =Y.^-^m\0?^{i^K)\QQ) • (6) 

pert QCD m'^"^^ port NRQCD 

n 

When quark and antiquark are in a particular angular momentum state J and color state 
1 or 8, the imaginary part of the left-hand side of Eq.(|6]) can be calculated with the covariant 



projection method [lO|: 



pert QCD J\ J '—^ ■' 

(7) 

where {0{'^^^^Lj))qcd equals to the corresponding NRQCD four-fermion operator expec- 
tation value at tree level, L is the orbital angular momentum, S is the total spin of the 
heavy-quark pair, and K is the degree of freedom of the initial state. For spin-triplet states 
with L=0, L = 1, and L = 2, the relations between M.{QQ)\^f'j and the full QCD Feynman 
amplitude A4 are 

MiiQQfsl ^ LH) = e^^^Tip^'^^WM] |,=o, (8a) 
MiQQfl, ^ LH) = eg/] ATr[C[«]n''A^]|,=o, (8b) 

MiQQfh^^ ^ LH) = ^4^/p'd^Tr[CWn''M]|,=o, (8c) 
where e^, e^-l and e^i are the polarization tensors for L = S, P, D states with total angular 



-pi ap 

momentum 1, Pj, Dj, respectively. For spin-trip 
heavy-quark pair accurate to all orders of f ^ is 



et states, the spin projector of the incoming 
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where is the four momentum of the heavy meson and P^ = AE"^, and 2q^ is the relative 
momentum between the quark and antiquark. The color projectors are cm = ^ and 

In the Fock space, the ipi^Dj) states are represented by 

\'Dj) = 0(1)|QQ(3dW)) + 0{v)\QQCPf^)) + 0{v')\QQCS^j^'^)) + ■■-. (10) 



Here, the probability of P-wave and 5'-wave states are suppressed by f ^ and relative to 
the D-wave state respectively, but their operators scale as and f""^ relative to Oi{^Dj). 
The relativistic effect and other Fock state contributions are suppressed at least by w^. Thus 
at leading order in v"^ the NRQCD formula for '^Dj decay into LH is 

r ("Dj ^ LH) = 

where the four-fermion operators are defined^ as 





(12a) 




(12b) 


Osi'Po) = I^Hy"^ ■ ^)t^xx\y'^ ■ 


(12c) 




(12d) 




(12e) 


OiCD,) = Y^fK^X X^K'^, 


(12f) 


0,('D2) = ^i^^K^'x X^K'^^, 


(12g) 




(12h) 


are = ct^ S'^ , K^^ = e'^^(T^S^^ + e^^^a^S'^, K^^'' = 





fjkgij _ 2^i^s^kgii ^ ^kigji ^ ^iJski)^ where S'^ = {^ fC^'"^^ - l^^6'^). 

For some processes, we need to calculate the NLO QCD corrections. To handle the 
ultraviolet (UV) and infrared (IR) divergences in the dimensional regularization scheme, 
one should extend the projection method into D = 4 — 2e dimensions. The definitions of 7 
matrixes in D dimensions can be found in quantum field theory books. And the sums over 
polarization tensors ejf^, e^ap^ and e^^^^ in D dimension are 

E^^4'^* = npP'. (13a) 

Jz 

^ The normalizations of the color singlet four-fermion operators agree with those in Ref.[lo| 
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E^S4v = 7ATn.pn.v', (13b) 



"ap a' p' 



E ^S^iv = J(n-'npp' + n^na'p) - Tr^n^.n^,,', (i3d) 



2 \ UU pp • ULp U p, JJ — \ 



^Sp^a'/jv - 2(D - 2)(z:' + 1) (ii"pii"V'n/3/3' + n/spn/j/p/iiaa/ + napn^s/p/n^'/? + ^fj^^a'p'^af^' 

J z 

2 4 

— -J^ -^J^ap^aip,'^Ppi + n^pHa'/j/nap' + nQ/p/IIa/Sn/j/p + n^/p/IIa^na/p) + -jj^ -^^i^a^Tlai pi\ippi\ 

(13e) 

X^^o/3p^"'/3'p' ~ g(2nQ,a/n/3^/npp/ + 2na^/nQ//3npp/ — nca/n/jp/iip^/ — no^/n/jp/iipct' — Hapi^^^'^pa' 

J z 

— nQp'n^Q'np/3/) + -^{—^a^^a'^'^pp' + Snap'IIa'/J'n^p + Ulap^lipi^a'pi + sna/siia/pn^/p' 

+ 2nQ,/3n^'pn„/p/ — nQ,/3'nQ,/p/n^p — ncpiia'/jn/j/p/ — Haa'^^s^^' pi — ^ap^^^'^ai pi)-, 

(ISf) 

E^^a/3p^a'/3*'p' ^ g (Haa' 1^/3/3' Hpp/ + nc^a/n^p'IIp^g/ + Yici^iYi^a'^pp' + na^/n^p/EpQ,/ + Yicip'^^P^'^pa' 

J Z 

+ nQp/n^^/iip^') — -^^^j^-—^{jiapYipaiTipi pi + na/3np^/no/p/ + ria/jiipp/nct'^' + Hap^^o'^p' p< 

+ napn^^/Ha'p' + Tiap^Ppi^a'P' + n^gpn^Q/n^/p/ + n^gpn^/S'IIa/p/ + Xi^f^ap'^a'^')- (13g) 

And the degrees of freedom are D — \ for the S-wave state; 1, (-^-iK-P-^) ^ ^^j, 
J=0, 1, 2 P-wave states; and D - 1, (^-3)(^-i)(^+i) ^ (z.-2)(d i)(d+3) j = 2, 3 D-wave 



states. The rather trivial extensions for L = and 1 cases were given in Ref. 10|. Here 
three principles are adopted to construct the nontrivial results for the D-wave case. First, 
the symmetry of the three indexes should be kept; second, the inner products between one 
tensor and the other two are zero; third, the completeness condition should be satisfied: 

E(l) (1)* (2) (2)* ^ (3) (3)* _ 

^aPp^a'P'p' + 2-^ ^aPp^a'13'p' + 2-^ ^ajSp^a' fS' p' ~ 

J z J z J z 

(-(naa'II^a^/ + liap'^a'p) — ^ _ ^ ap^a' ^'^^pp' ■ (14) 
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We calculate the degrees of freedom of each j state with group theory and get Eq. fll3ep . 
Then we derive out Eq.pSgP with the help of the first two principles. In the end 



Eq. (jm =Eq. dUD-Eq. (^^-Eq. ([Bg]) . 



To get the NLO NRQCD results, the operator mixing equation f[T^ between P-wave 
and Z^-wave operators in momentum space should also be extended into D dimension for 
consistency: 

J2Cpj,DAiQQ>D,,\OCDj,)\{QQ)sn,) = {{QQ>Pj\0('Pj)\{QQ)sp,)q- q', (15) 
where Cpj^D^, are the generalized Clebsch-Gordan coefficients: 

^^-^■^ = ,(^.) JI'/J'P'P.) ' ^^^^ 

where the repeated indexes mean being summed in D dimension. 

After resolving the above problems, one can do calculations from both full QCD and 
NRQCD straightforwardly. 

III. FULL QCD CALCULATION 

In Sec. II, it has been explained that at leading order of f ^ the LH decays of ^Dj contain 
the subprocesses of QQ^I^g^, QQfpj{J = 0,1,2), and QOtl^jiJ = 1,2,3) annihilating into 
gluons or light quarks. When doing the calculation with full QCD theory, for simplicity, 
only the explicit imaginary part of A is given here, 

2W((QQ)Kl^^ ^ (QQ)Kli,) =1? IY. \M{{QMAl, ^ LH)|2rf$. (17) 

pert QCD A J '—^ ■' 

In the following subsections, the contributions at order, and the corresponding real and 
virtual corrections will be given in Secs.IIIA, IIIB and IIIC respectively. As for the processes 
whose tree level diagrams are already at a\ order, their results will also be given in Sec. III. 

A. LO Results 

— [Hi ^ f8l 

There are three subprocesses {QQyA^^ — )■ gg, {QQyip^ ^ — >■ gg with nonvanishing imaginary 
parts at 0{a1). The Feynman diagrams are shown in Fig. [1]. And the results in D dimension 
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FIG. 1: Feynman diagrams ioi{QQ)\,„ — )• qq, {QQ) 



Si 



"■PO.2 



99 




— roi 

FIG. 2: Feynman diagrams for {QQys^^ — )• qqg 
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(2W(^SfO) 



3c[8hABorn_ 8iVXvrV^$(2) 



3-2e 



Q 



3-2e 



(2Im^(^Pf )) 



[8K,Born_ 4e^ " 13e + 6 32i?paf7rV^$(2) 



where $ 



(2) 



Stt r(2-2€) ^mi 



(3-2e)(5-2e) 
is the two-body phase space in D dimension. 



[18a) 
;i8b) 
(18c) 



B. Real Corrections 

Besides the real corrections (RC), other processes with three bodies in final states include 
iQQfs'! ^ 3(7, {QQfl^ ^ ^9,{QQfl^ ^ qqg and {QQ^l^ 3^?. And the typical 
Feynman diagrams for qqg and 3g are shown in Fig. [2] and Fig. [3], respectively.^ 

We calculate the Feynman amplitudes and do the phase space integrations both in D 
dimension, and check the results with eikonal approximation relations or Altarelli-Parisi 



^ In Feynman gauge the ghost diagrams are also needed, where three-gluon or four-giuon vertex will appear. 
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FIG. 3: Feynman diagrams for {QQ)3j[i,s] — >• 3g. Diagrams with different positions of ki are 
neglected. In Feynman gauge the ghost diagrams should also be contained. 



splitting functions when one gluon is soft or two light partons (gluon or light quark) are 
coUinear. Three variables Xi = — are introduced to describe the (QQlsWi.si ki + k2 + 
process. For energy conservation Xi + X2 + = 2. The three-body phase space in D 
dimension becomes 

4:771% TT £ 

(i$(3) = $(2) . N2rn - [(l-a:i)(l-a:2)(l-a^3)] 5{2-xi-X2-x^)dxidx2dx'i. (19) 

In the phase space, the Xj = region is the soft region of particle with momentum fcj, and 
the Xi = 1 region is the collinear region of the other two particles. And the inner products 
of those four momenta are P ■ ki = 2mQXi, ki ■ kj = 2mQ{xi + Xj — 1). 

When L = S and P, the real corrections contributing to 2Im^(^Lj'^')^*^ are 



1 



^ iK.,M) lT.\Mi{QQ)l!fl ^ qqg)\'d^i3), (20) 



where Kd = 1, Kcs = — 1, and Kj is the polarization number of angular momentum J 
state. In the D-wave case, because of parity conservation, only the 3g process is left. 
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^- {QQ)3T{m 999 
'-'J 



There are 18 diagrams in the QQ — > ggg process. Because of J^*" conservation, only class 
[gl] is needed in ^S'f ' — > ggg and ^Dj' — > ggg processes. In the ^Pj^' — > ggg processes, the 
diagrams are those in class [gl], class [g2] and the corresponding ghost ones in class [Gl]. 
And all the diagrams should be calculated in the ^S'f — )■ ggg case. We now proceed to 
show how to get the physical results at a differential cross section level when ghost diagrams 
are contained. In the Feynman gauge, when the three or four gluon vertex appears, the 
nonphysical (NP) degrees of freedom are removed by ghost diagrams. Label the square of 
the amplitudes of the gluon diagrams in the Feynman gauge with the subscript NP, and 
express the ghost result as 

Wg\U, = E ^ 9GHfi^,)\\ (21) 

where G^. and G\^. are for ghost and antighost with momentum and /cj, respectively. To 
cancel the NP part of each gluon, all possibilities of indexes should be summed over. 
Then proper physical results at the differential level are 

E miOQtt] ^ 999)? = E \^^PmQtt] ^ 999)? - E (22) 

Here we omit the details and only give the real corrections below 

E \Mm)fl ^ 999)?d^^s, = '^'TnJf^"^ (23b) 



3!(iV2 - l)(3-2e) y ''^^^'^1 108m, 



3,(^2,1) J E IMiQQfk ^ 999)\'d^i3) = 

^'(^ + I + ^I^)(2W(3pj«i))Born, (23c) 



-32i) / E \MiQQ)fk ^ 999)\'d<^>i 



_ CABpi-lSSir^ + 1369)af 

3\{N^ - 1){3 - 2e){2 - 2e) J ''"vv^^y.Pi ' y^^yi -^(3) - , 

(23d) 
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3!(3 - 2e) 

32 ^ ^ ^ s 4./ 608 -7744 + 16057r\ , 
^^B,^.,F,.ay(-— + ^^^^) (23f) 



3!(5-2e)(3-2e)(l -2e) 



32 „ „ 3 4 8 , -23024 + 21257^^ 



3^^i3AF,.„i;,-(-- + ) (23g) 

3!(7-2.)(3-2.)(2-2.) J T.\mQQto, ^ = 

B.^.FM,"(-- + + (23h) 



3m^ ^ V . g3QQ 



where the are correct at = (^) r(l + e) ,Ca = 3, and Bp = 5/12 are color 

[1] 



factors. The expressions of ^ \M.{{QQ)tl) — > 999)1'^ are too comphcated to be given here. 



2- (<3Q)3r[i.8l QQ9 



At C(q;s) there are four subprocesses: {QQ)^^^ — )■ qqg, which include all the Feynman 
diagrams in Fig. [2], and {QQ)3pi8] — >■ qqg, in which there are only two diagrams q[2] and 
^[5]. As in the 3g processes both the Feynman amplitudes and the phase space integrals are 
calculated in D dimension directly and the results are 

(Ar.-iK3-2.) / E IMiiQQfl ^ qq9)\'d^i3) = 

^ {Cf{^ + Ye + t) + + ^)) (2W(35pl))Born (24a) 

PfTT) / E |A^((W)6o ^ QQ9)\'d^i3) = 

(^/^(-^))(2We^(l'^))^°'"'^ - ^5F«.(3(2W(3M«l))Born|. + ^^2^ (24b) 

(]vFi) / E |A^((W)il ^ qq9)\'d<^>i3) = 
(^/^(-Ji))(2Ini^('A^''))^°^° - 94S^«,(3(2W(35fl))Born|k + l^al) (24c) 
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XNi^IE\Mi{QQ)fl^^qqg)\'d^^s) 



C. Virtual Corrections 

The virtual corrections are performed with a renormahzed Lagrangian, in which the 



renormahzation constants of the QCD gauge couphng constant gg = y/Anas, heavy-quark 
niQ, heavy-quark field ipg, light quark field ipg, and gluon field are defined as 



g^s=Zggs, m% = Zm^mQ, ^1 = ^q^, = v^V", AJ = v^^^, (25) 

where the superscript labels bare quantities, and Zi = 1 + Here a mixing renormal- 
ization scheme jl^ is adopted. The quark mass mg, heavy-quark field i/jq, light quark field 
ipn, and gluon field are defined in the on-shell condition, while gg is in the minimal- 



subtraction(M5') scheme. Then in this mixing scheme, these renormahzation constants are 



5Zgf = -^£^^±- + A + 3 ln(4) + 4 + 0{e)) (26a) 
^ 4:71 e„„ e 



^^os ^ -^^( A - ± + 0{e)) (26b) 



IR 



6ZOS = ^(/3o - Ca){— - — + 0{e)) (26c) 

47r err,r £ , „ 



^os SCpasF,, 1 ,,,,,, 4 



5^™, = + 1^(4) + o + 0{e)) (26d) 



SZr-^i^-H^) + Oi.)) (26e) 



7MS 



where /3o = ^^y^ — and Nj is the number of light flavor quarks. The representative 
virtual correction Feynman diagrams for the Born processes in Sec.IIIB are shown in Figs. 4 
and 5, without external leg correction diagrams in this scheme. The UV divergences in self- 
energy and triangle diagrams will be canceled by the counterterm diagrams correspondingly. 
To regularize the Coulomb ^ poles in the virtual processes, the loop integrals are done first 

— [81 

before setting the relative momentum g = 0. Also in {QQ)ap^^ — ?■ gg processes, we integrate 
the loop momentum first then compute the first derivative of the Feynman amplitudes with 
respect to q"'. When q ^ 0, the momenta of Q and Q are Pq = ^ + q, Pq = ^ — q, where 
P is the meson total momentum, and the momenta of the two massless final state particles 
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FIG. 4: One-loop Feynman diagram for {QQ):^ pis] — >• gg 

"0,2 

are labeled with ki and k2- Then the scalar functions of the loop integrals can be expressed 
in the Mandelstam variables: 

. = (PQ + PQ)^ t=iPQ-k,f, U={PQ-k2)\ (27) 

where s = 4mQ/(l — f^). We do the calculations diagram by diagram and summarize the 
results in the following form: 

2lmAeL^]''Y'' = 2 W(3l5'^])^°-^ y D,, (28) 

k 

where for each process are listed in Table[I], [II], [III]. The contributions of the counter- 
term diagrams are put together with the corresponding self-energy and vertex diagrams to 
show explicit cancelation of the UV divergences. The IR divergences left will be canceled 
by the RC. 
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FIG. 5: One-loop Feynman diagram for {QQ)sol»] ~^ QQ 



TABLE I: Virtual corrections to {QQ)3„[s] — )• qq- 



Diag. 


Dk 


61 + 64 + 66 


Ca-2Cf I Ca-3Cf-I3o , Ca{9-27t^)+4Cp{-6+7t^) /3o . / fi^ \ 




62 + 63 + 65 


i2CF-CAW Po , CA(7+21n(2))-12Cp , ^„^ M^" ^ 
4. 2e,^ 1 3 1 2 


67 


LB. 


68 


(2C^-Ca)(^-4) 

IR 


69 





TABLE II: Virtual corrections to (QQ)3p[8] gg- 



Diag. 


Dk 


al + a2 


Cf I Cj.(5+321n(2)) 


a3 + a4 + a9 


-Ca , -9A)-Ca-18Cf , CF(-44-1281n{2)+37r2) , Ca(38+8 ln(2)-7r2) , « ^2 


3e2^^ ' 9e,^ I 18 1 36 1 '^uin^^^^j 


a5 


CaI \ 2 14 1 13 „2 1 14 
2 ise^ 1 9e,„ 27 + 18^ + 27 ^^y^>) 


a6 





a7 




a8 


Ca{^ + 'SH2)) 
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TABLE III: Virtual corrections to {QQ)sry[8] gg- 

9 





^ K 


al + a2 


Cf(^ + ^ + ^ln(2)) 

^ t V e '66 \ // 
in 


/j 3 _|_ /j4 _(_ 


-3Ca I -32/3o+15Ca-64Cf Cf(100+172 ln(2)+37r2) (75+112 ln(2)+187r2) ^ 1^/ \ 


16€^ ' 32e,„ 24 ' 192 ' ^'^Um^ 

J ff IK Q 


ah 


Ca( 13 15 17 1 37 „2 1 89l^/'o^^ 
2 I 8e^^p 16e^^ 288 + 48^ + 18 ^^y^l) 


06 





a7 


(^^-%)(^ + 2^-'-f + |-' + iM2)) 


08 


-C^(| + f ln(2)) 



D. Summation of Real and Virtual Corrections 



Collecting the RC with VC, we obtain the full QCD results at O^a^)- 



(2Im^(=^,sW)) 



40^3 (tt^ - 9) 



pert QCD 



81m? 



(29a) 



(2W(34«])) = ^(5a,(-657 + GTtt^) + Ar^(367r 

pert CJCD Q 

+Q;,(36(2Ci. - Ca)£ + 642 - 207V^ - 2777^ + 72 ln(2) + 36^o In ^gr))), (29b) 



{2lmA{^Pr)) 



5a? 



pert QCD 432™'^ 



^(2167r + a,(54(2C^ -Ca)- + 3032 + 2l7r2 



+8401n(2) + 216/3oln(^))) - 3^S^a,(3(2W(34«]))Born|k + g^^^^^, 



(29c) 



(2W(3pfl)) 



9m: 



_ 5a^(1369-1387r^) 
pert QCD Q 

4-S^a,(3(2W(35f ))Bo-|. + '^a% 



(29d) 



(2W(3pj«])) 



pert QCD 



2,^(1447r + Q;,(36(2Cf - C^)^ + 4187 - 2587r2 



+3361n(2) + 144/3oln(^))) - 5^S^«,(3(2W(34«l))Borng + ^^2)^ 



(29e) 



(2W(3dW)) 



32 



3,,4e/ 



608 -7744 + 160577^ 



pert QCD ^TTIq ^ ^ " sr \ -j^g^^ 



+ 



16200 



) (29f) 



(2W(3dW)) 



32 



pert QCD Bmg 



D ^ 17 3 4., 8 , -23024 + 21257r2 



15e 



1800 



) (29g) 
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(2W(^dW)) = -lZ^Bp<!>,F,7iaifi''^{-^ + "^""r^;, ) (29h) 

port QCD ottiq 15e 6300 

There are still infrared divergences and Coulomb singularities in some of the expressions 
above. As explained in Ref.j6|, the infrared divergence comes from the soft gluon emission 
of heavy quarks, and the Coulomb singularity reflects the behavior of heavy quarks in the 
potential region. In next section, both of them will be repeated precisely when doing the 
NLO corrections for NRQCD matrix elements in the corresponding regions. 

As mentioned above, the S- and P-wave subprocesses have been studied by Petrelli 



32 „ ^ ^ ... 32 -28656 + 26457r2. 



etaL[lO|]. In their paper, the soft and collinear singularities are separated with the help 
of eikonal approximation and Altarelli-Pasrisi splitting functions, then they calculate the 
finite part in 4 dimension. In this paper, we recalculate them in D dimension directly as a 
cross check, and get the same results. The D-wave subprocesses have also been considered 
in Refs.y, Q] but they did the calculations in 4 dimension, and regularized the infrared 
divergence with the binding energy. 

IV. NRQCD RESULT AND OPERATOR EVOLUTION EQUATIONS 

There are three typical energy scales in the heavy quarkonium system, related to the 
small parameter v. They are mg (the heavy-quark mass), mgf (the typical momentum of 
heavy quarks in heavy quarkonium), and rnqv'^ (the binding energy). Then, there are three 
dynamical regimes in the NRQCD effective theory, in which either the heavy-quark or the 
gluon is on mass shell, and they are 

soft regime : : A;o ~ \k\ ~ mqv, : T ~ |p| ~ rnqv 

potential regime : : ~ ^qv'^, \k\ ~ mqv, : T ~ mgv'^, \p\ ~ mqv (30) 

ultrasoft regime : : fco ~ ~ rnqv"^, 

where fcj, and Pi, are the momenta of the gluon field and heavy-quark field, respectively, 
and T = pq — mq = + 0{v'^). Because there are more than one regimes in the non- 
relativistic system, matching the production and annihilation of external heavy-quark and 
antiquark pairs at certain order in v can not be manifest, though the power counting rule, 
velocity scaling rule, of operators in NRQCD is simple. This problem has been addressed in 



several papers 18N25|. and the matching prescriptions based on dimensional regularization 



in NRQCD were also clarified. Furthermore, the potential NRQCD (pNRQCD) effective 
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soft 



{T,p) 



iT,P) 

potential : Q,, : » 



k J 



all 



ultrasoft : Au : 'WWW = 



FIG. 6: NRQCD Feynman rules for heavy-quark and gloun propagators in different regimes. 



(T., 



(a) 



(6) 



FIG. 7: NRQCD Feynman rules for heavy-quark gluon vertex. 



theory 



Ref. 



was proposed by introducing the potential to manage the nonperturbative effect in 



26| . The NRQCD Feynman rules for propagators in the Coulomb gauge|25|], are shown 
in Fig. [6], where = 5ij — ^p-. At leading order in f^, the interaction term of the NRQCD 
Lagrangian , where the antiheavy terms are neglected, is igsi^^AQ + And it 

turns to be igsil^'^{AQ + ^^)ip^ for V • A = in the Coulomb gauge. The Feynman rules 
for vertex can be read directly and are listed in Fig. [7]'^. The Feynman rules for anti-heavy 
quark can be gotten by charge-conjugation symmetry. 

Since the short- distance coefficients are obtained by matching full QCD results with 
NRQCD results, we only need to calculate the real parts of the matrix elements. Figure 
8 gives the LO Feynman diagram. At NLO in a^, when the inner gluon line joints two 
incoming or outgoing quark lines, a nonvanishing real part only appears in the potential 
region. When the inner gluon line connects with one incoming quark line and one outgoing 
quark line, the power counting rules 25| tell us that the soft region will provide the leading 
order contribution in v. The external self-energy diagrams are dropped to be in accordance 



The Feynman rules are the same for the corresponding interaction terms in different regimes though their 
power counting may be not. 
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FIG. 8: NRQCD Feynman Diagram for LO matrix elements 



(T,p) (T'.p') (T.pj {T\P) 




FIG. 9: NRQCD Feynman Diagrams for NLO matrix elements 



with the renormahzation scheme in full QCD calculation. Then we only need to calculate, 
two class, six NLO Feynman diagrams shown in Fig. [9]. 

For convenience, we present the detailed NLO corrections of the P-wave octet matrix 
elements, (C(^Pj^')), which are more representative than the S'-wave ones. The LO result 
{0{^pf^))Bom is trivial. Using the Feynman rules for propagators in the soft regime and 
vertices, the loop integral of diagram (a) reads 

tg', f d^k p-p'-(p-k)(p'.k)/k2 1 1 ^^^^ 



rriQ J {2n)^ fcg — + ie ko + ie ko + ie 

After performing the contour integrating of ko = \k\ — ie, 

^ _ gl f d^-'k p.p--(p.k)(p--k)/k^ 

which is both infrared and ultraviolet divergent. In the dimension regularization scheme the 



result is 



la = - -)P ■ P' (33) 
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The integrals of (b-d) in Fig. [9] could be calculated in the same way: 

a, . 1 1 



I b-d 



Svrmi, euv e 



-)P ■ P' (34) 



The loop integral of diagram (e) in potential regime could be written down similarly: 

= 1 jl^v^JTi;^i^TZ¥^l-^^^^. (35) 

where T = When ko is integrated out: 

2 f d^-'^ 1 1 
= '^"^^ J (2;0^k^ k^ + 2p.k-ze 

This could be done in D dimension directly. Regularize the Coulomb singularity by intro- 
ducing V = — , then at f order we have 

mo ' 



. a--(--ln(^)-7i.)) (37) 

The integral of diagram (f) has the same real part but with a plus sign before the virtual 
part. 

The color structures for diagrams (a), (c) and (b), (d) and (e), (f) are obtained by the 
color decomposition listed in the first, second, and third line below respectively: 

V2T'^T^ ® V2T''T^ = Ci.^ ® ^ + ^y2T^ ® V2T'', (38) 
rpb^T^T^ ® V2T^ = {Cf - \Ca)^T" ® V2T^ 
Combining the integrals with the according color factors and summing them over, we 
obtain the NLO NRQCD corrections for the P-wave octet operator matrix elements, which 
are UV divergent and need to be renormalized: 

{0'^{^Pf))NLo = {{l + ^{CF-\CA))V2T'^®V2T-+ 
4". _ 1 ^ 1 + BfV2T^ ® y2r^]p ■ v'}{Oi'Pj)) 



where 0{^Pf) = d{^Pj)V2T'' O V2T^, Bf = and the superscript "0" means the 

matrix elements of the bare operators. As expected, at NLO the P-wave octet operators are 
mixed with the D-wave singlet and octet ones, and with the help of Eq. ([T5l) . they could be 
reexpressed as 

m'P?))NLO = (1 + ^{Cf - '^CA)){OrPP)Uorn 
+ '-^i7^-'JiCF{OrD^i^)horn + BF{OrD^!}))Born), 
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where Cj^ji are defined in EqlTGland for J' = 1, 2, 3, C*o,j' = -^^^^i^t^5 0, 0; Cij/ = j^j^, |, 0; 
C2,j'= 4(^~^i)2 ; J, 1- The P-wave operators are mixed with the D-wave operators at NLO in 
as- But the NLO NRQCD corrections of D-wave operators are related to the relativistic 
corrections of P-wave operators. Then at leading-order of f^, the renormalization transfor- 



mations of those operators in MS scheme are 

OrP?) \ ^ I 1 C^;,(^ + ln4vr-7i.)\ / O^H^Pf) 
0\^D^]}) / V '''' 1 ) \ 0^{^Df, 



(41) 



where Cjj, = Cj^j'^^^. When the operators are mixed with each other, the renormal- 
ization constants Z' s are not numbers but matrices. The D-wave octet operators are also 
dropped, for they do not appear in full QCD calculations. The matrix elements of the 
renormalized operators O^i^Pf^) at NLO are now UV finite, but there are still \ poles 
in D-wave terms, Eq.( H2|) . and Coulomb singularities, which will absorb the infrared and 
Coulomb divergences in full theory: 

{OyPP))NLO = (1 + ^{Cr - lCA)){0^fPP))Born 

where /ia is the renormalization scale. The matrix elements of the S'-wave octet operator at 
NLO could be computed in the same way: 

^ m'Sf^))r,Lo = (1 + ^{Cf - iC^))(0^(^5p^))Bor„ 

The matrix elements of the S-wave singlet operator and D-wsve singlet operators at NLO 
do not need to be calculated , for their LO short- distance coefficients are already at 0{al). 

Multiply the matrix elements with the short-distance coefficients, we obtain the NRQCD 
result at NLO in ct^: 

{2ImAeS?)UnQCD = ^^"^i'^^'' '))Ln (44a) 
{2ImAesP)hnQcn = ^-^^(1 + ^(C, - ic^))(0(^M^l)i_ (44b) 



r97-^^r3p[8Kx _ ^ 2/m/(3p^') g.TT 1 
[^AimAy Pq ))nrqcd - { ^ (i + ^l*-^^ ~ 2 

4a,P^2/m/(35fl),4V ,,^,s ^[su.r 



^ .r(l + .)}(0(^Pr))g.„ (44c) 
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- ^^^^^(^)'r(l + ^)}(OePf' ))g,„ (44d) 

-f8K 



t!:^^^'"»^'' (^)T(l + (44e) 



37rml e //i 



mi 

E '°w"'""f''''' (^)T(l + .)}(OeM"))f«„ (44£) 



m, 

- E + (44g) 

- E '°w"'""T''''' (^)-^" ^ ^))(oe4"))L„ (44h) 

Finally, matching the NRQCD results with full QCD results, we get infrared safe short- 
distance coefficients at 0{a^): 

2Im/('S|'l) = 12^1^, (46a) 



2Im/(''sf ') = :^(36A'/7r + a,(5(-657 + 677r') 

iUo 



+ iV/(642 - 20^"/ - 277r^ + 72Zn2) + 72/3oNfln-^)), (45b) 

2mQ 

r 2 

2Im/(^P(j^l) = ^^(6487r + a, (9096 - 464A^/ 
1296 

+ 637r^ + 2520/n2 + 1296/3o/n-^ + 967V//n^^)), (45c) 

2mQ iiA 
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2 

2Im/(^Pi^^) = ^(4327r + a, (12561 - 464 A^y 



- 774:71^ + 1008/n2 + 864:f3oln-^ + 2A0Nfln'^^)), (45e) 

2mQ /XA 

(32l7r2 - 8032 - 29184:ln^)a^^ 
2lmfCD?) = (45f) 

(4257r2 - 4816 - 384/n:5^)a3 
2lmfCD?) = — (45g) 

(5297r2 - 8688 - 5376ln^)a^, 
2Im/(^Z^W) = ^ (45h) 

The P- and D-wave short-distance coefficients are /ia dependent, and their /ia dependence 
can be canceled by the renormahzed operator /^a dependence. The /ia dependence of the 
renormahzed operators could be derived out by finding the derivative of both sides of Eq. ( l4Tll 
with respect to /ia: 



J' 



(46b) 



dO^^^^ _ y _^l^o^{'Pf) (46c) 



d In /iA J 2Nc 37rm| 

Remember, for bare quantities, — — - = 0. For a phenomenological reason, we also give 
the /iA dependence of C(^S'f^), though we do not calculate its NLO NRQCD corrections. By 
solving the differential equations, all S- and P-wave operators' expectation values in \Hji) 
states are related to that of the D-wave singlet operators: 



in^^{H,.\OrD'j))\Hj,), (47a) 
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2 3m^/3o as(/iA) 

In pNRQCD, the S'-wave color-octet matrix elements for P-wave heavy quarkonium de- 



27, 



And the relations 



29| . which shows that the 



cays are also estimated through operator evolution equation 
between their results and ours are discussed in our previous work 
two methods are consistent with each other. 

V. NUMERICAL RESULT AND DISCUSSION 
A. ■^Dj Decay into LH 

For heavy-quark spin-symmetry, the long-distance matrix elements of D-wave four- 
fermion operators are equal to each other for different J, and relate to the second derivative 
of wave functions at the origin: 

{H,\0{^Df)m = {Hs\0{'D^i^)\Hs) = Hnml 

The matrix elements of the P-wave octet operators and the S'-wave singlet as well as octet 
operators in the corresponding J' states could be estimated through the resolution of op- 
erator evolution equations, Eq.( H7|) . When /ia^ and /ia are separated widely enough, the 
evaluation terms will be much more important than the boundary terms labeled with /iAo- 
Here we set /xaq = "^Q"^; where = 0.3 for charmonium and = 0.1 for bottomonium, 
since the NRQCD perturbative calculations could only hold down to scale of order mgf : 



{HJ,\0\'pf){^,^)\HJ,) = Cj/^\n^^Hnm'Q, (49a) 
{HAO\^Sf){^^^)\HJ,) = ^(4ln^)^i/z.m^, (49b) 
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{Hj,\0^eS?){^.,)\Hj,) = £^{^Jn-^rHnml, (49c) 

We also assume /xa = /i, for the factorization scale /xa in NRQCD also acts as the renormal- 
ization scale in operator renormalization. In the end, we come to the overall expressions for 
the LH decay widths of ^Dj{J = 1, 2, 3) states to NLO in as at leading order of v : 

T(^Di LH) = (-5.0«^(0.167 + In -^) + 

2mQ 

«2(15.7 + a,(88.5 - 4.33iV^) + a,(10.0/3o - 1.32iV^) In ^) In ^ 



/3o 

1.32al{1.57Nf - 0.278«,(iV/ - 21.4) (iVj + 0.093) + 1.0iV/a,/3o In ^) In 



2 as 



A 



2 )Hd. (50a) 







TCD2 LH) = (-0.59a^(1.62 + In ;r^) + 
a2(1.87 + a,(8.14 - 1.95iV/) + a,(1.19/3o - 1.32Ar^) In ^) In f 



/So 

1.32a2(i.57Ar^ _0.278a,(iV/ -21.4) (A^j + 0.093) + 1.0A^/a,/3o In 2^) In^ f 



2 



(50b) 



r(3D3 ^ LH) = (-2.37a^(0.645 + In -^) + 

2mQ 

a2(7.45 + a,(30.8 - 2.55iV^) + a,(4.74/3o - 1.32iV^) In ^) In ^ 

Jo ^ 
1.32a2(i.57iV; - 0.278a,(iVj - 21.4) (A^j + 0.093) + 1.0A^/a,/3o In ^) In^ 

2 ^)-?^j:> (50c) 



Ao 

where = "^(/^Ao)- 

1. D-wave Charmonium D j) LH Decay 



Making a choice of = l.SGeV, Kqcd = 390MeV, Hdi = = 0.786 x lO-^GeVjaOl 

and Nf = 3 for charmonia, we obtain at /i = 2mc: 

T{tlj{l^Dj) ^ LH) = (435, 50, 172)keV for J = (1, 2, 3). (51) 
When fi = rric and the other parameters are fixed, the results turn to be: 

T{iIj{1^Dj) LH) = (683, 42, 223)keV for J = (1, 2, 3). (52) 
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FIG. 10: Renormalization scale dependence of the decay widths of charmonium states I'^Dj to LH 
at order. 



And the yU dependence of the decay widths at 0{al) is shown in Fig. [10] 

In the potential model, iIj{1^Dj) can only decay to 3g at order. The infrared diver- 
gences are regularized by e^e the binding energy of the bound states. Accurate to ete order, 
potential model results|3|] are: 

TcCDj LH) = (160, 12, 68)keV, for J = (1, 2, 3). (53) 

If we reset their parameters the same as ours with as = as{2mc), M = 2mc = 3.0GeV, 
= O.OlSGeV^, the potential model predictions become: 

TcCDj LH) = (240, 18, 102)keV, for J = (1, 2, 3). (54) 



It could be found that in the cc system the NRQCD predictions are about 2 ~ 3 times larger 
than potential model results. In leading logarithm approximations [4], the ratios of the LH 



decay widths for J = 1,2,3 are T{^Di) : T{^D2) : T(^D3 
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1:4. Including the non- 



negligible corrections to the leading logarithmic terms [3|, the ratios turn to be: 40 : 3 : 17. 
And the relative ratios predicted by NRQCD at /i = 2mc = 3.0GeV and /i = rric = 1.5GeV 
are 43 : 5 : 17 and 34 : 2 : 11, respectively. 

Much works has been done to predict the mass spectrum of ipll^Dj); some of the numer- 

32I, and some theoretical work reviews may be found 



ical results are collected in Refs. 



31 
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in Ref. 33| and references therein. All the predictions indicated that the masses of ip{l^Dj) 
are all larger than the threshold of DD (about 3730MeV) , and the center of gravity of 
ID states calculated in the Cornell potential 34 1 is 3815Mevj35|. For its decay to open 



charm to be kinematically allowed, tpll^Di) should not be a narrow state. I t is be 
that '?/'(3770) is primarily a l^Di state with a small admixture of the 2^5"! state 
the latest experimental average of its width is T{ip{3770)) = 27.3 ± l.OMeV 



ieved 
36|,|37|, and 
. But there 



is a long-standing puzzle in its non-DD decay that the '?/'(3770) decay is not saturated by 



the DD decay 



paper 



39| . A detailed discussion about this problem could be found in our previous 



40|, and in this paper we will briefly review it in Sec.VB. 



The remaining J = 2 and J = 3 states are both expected to be narrow with different 
reasons. iIj{1^D2) is presumed to lie between the DD and DD* thresholds 4l|] and is forbid- 
den by parity to decay into two pseudoscalar D mesons. While narrowness ofj/'QADa) in 



contrast is due to suppression by the DD F-wave angular momentum barrier |31l. I4l| . The 
principal decay modes of ip{^^D2) are radiative transition {ip{l^D2) — )■ 7Xci5 7Xc2), hadronic 
transition {iIj{1^D2) — > J/ipmi), and LH decay. To ip^l^D^), these decay modes are also 



considerable since r{ilj{l^D^) DD) is predicted to be only about O.SMeV 42|. when its 



mass is 3868MeV. And the decay widths predicted by the model 4l|, |42| including the 
influence of open-charm channels are r{ilj{l^D2) — ?■ 7Xci) = 212keV, r{ilj{l^Dj) — )■ 7Xc2) = 
(45, 286)keV, for J = (2, 3) and T{tp{l^D3) DD) = 0.82MeV at m^^i3n2) = 3831MeV and 
m^(i3B3) = 3868MeV. They also estimated T{^p{l^Dj) J/^pmr) = 68 ± 15keV. Using the 
numerical values in Eq.f lSTl) . we then roughly predict that the branching ratios for the LH 
decay of ip{l^Dj) are: 



Bt{iIj{1^Dj) LH) = 13.3%, 13.3%, for J = (2, 3) 



(55) 



2. D-wave Bottomonium T(n^L'j) LH Decay 

Unlike charmonium, T{n^Dj) (for n=l,2) are predicted to lie below the BB flavor thresh- 
old, and expected to be quite narrow, where n is the level number. Some predictions of 
T(l^Dj) and T{2^Dj) masses are reviewed in Ref. I43I . Taking = 4.6GeV, Aq^ = 
340MeV, Nf = 4, Hm = ^^^^ = O-^Ol x IQ-^GeV for ID states, and Hd2 = ^gfffj' = 
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FIG. 11: Renormalization or factorization scale dependence of r(T(l^Dj) — )• LH) at order. 
The solid, dotted and dashed are for T(l'^Di), T(1^D2)) and T(1'^D3), respectively. 

I — I 

0.750 X 10-^GeV for 2D states^], at /i = 2mfe, we find 

T{T{l^Dj) LH) = (6.91, 0.75, 2.75)keV for J = (1, 2, 3) (56a) 

T{T{2^Dj) LH) = (12.9, 1.40, 5.14)keV for J = (1, 2, 3) (56b) 
When /i = mfe and the other parameters are unchanged, our predictions turn to be 

r(T(l=^Dj) ^ LH) = (7.99, 0.60, 2.85)keV for J = (1, 2, 3) (57a) 

T{T{2^Dj) ^ LH) = (14.9, 1.21, 5.33)keV for J = (1, 2, 3) (57b) 

The II dependence curves of T{l^Dj) and T(2'^Dj) LH decay widths are similar, so only the 
n = 1 results are shown at 0{al) in Fig. [11] as an illustration. 

In the potential model, Belanger and Moxhay[4] found, for J = (1,2,3), the leading 
logarithmic results are r{T{l^Dj) — > ggg) = (2.2,0.26, l.l)keV, and a good approximation 
to the exact phase space integration given by Bergstrom and Ernstrom[3] brings a factor 
of 2 ~ 3 enhancement, and their results are F(T(1'^ i^j) — )■ ggg) = (6.3, 0.51, 2.7)keV, for 
J = (1,2,3). If we normalize them using our inputs at /x = 2mb and setting M = 2mh, 
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potential model estimations are then 



r(T(l^Dj) ^ LH) = (5.4, 0.51, 2.3)keV for J = (1, 2, 3), (58) 

which, to some extent, are in agreement with our NRQCD numerical predictions with /i = 
2mb- In the T{1D) case, it can be easily found out that the potential model results are 
dominated by the logarithmic terms. And numerically, the NRQCD results are mainly 
from the P— wave color-octet subprocess contributions. If we relate the logarithmic term 
ln(l/e) in Eqs. (20-22) of Ref. 3] to the evolution term In ^^^^^^ in this paper by setting 
ln(l/e) = In ^^jj^, we find the logarithmic terms as well as the vr^ terms in the potential 
model results can be exactly reproduced within the NRQCD approach. This then provides 
an alternative way to relate the value of In ^(^^°)^ to the potential model estimation. Using 



3], mft = 4.6GeV, as = 0.18, and Nj = 4, we get 



the inputs (r) = 2.5GeV ^ given in Ref. 
In "a^^^A) = 0.58, which is consistent with the value we obtained by choosing /iA = 2mi, and 
A^Ao = mbVh 



In Ref. 



44j . the branching ratios of some decay modes of T{l^Dj) are summarized in 



Table IX, where r(T(l^D^ ) — )■ e+e ) was calculated in Ref. 45| and r{T{l^Dj) — t- tttt) was 



obtained by Moxhay 461]. Since the LH decay widths of T{l^Dj) are now calculated in the 
framework of NRQCD, we update the theoretical predictions for these ratios in Table IV, 
where the numerical results in Eq.f l56ap are taken as estimations for LH decay widths of 
T{l^Dj). In 2004, the CLEO Collaboration observed T{1D) in the four-photon cascade 
process T iXb{2P), Xb{2P) ^ ir{lD), T{1D) ^Xb{lP), Xb{lP) ^ 7T(15), followed 
by T{1S) 1+1-, and the branching ratio is -6(7777/+/- )t(id) = 2.5 ± 0.5 ± 0.5 ■ lO-^ll?!. 
The signals are interpreted as predominantly coming from the production of T(1^D2)- Small 
contributions of T(l^Di) and Tll^D^) can not be ruled out. In the near future, with more 
accumulated data, all the spin-triplet T(l^Dj) states may be identified. Unfortunately, 
the D-wave bottomonium LH decays could not provide a good probe to find out whether 
NRQCD is prior to the potential model to describe the bottomonium system, for the differ- 
ence between the two theoretical predictions is small, unless a very precise measurement is 
made. 

For the n=2 states, no experimental evidence has been observed until now. To make a 
theoretical comparison for r(T(2^Dj) — t- LH), the numerical potential model predictions 
are needed. 
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TABLE IV: Summary of the partial widths and branching ratios(i?) for spin-triplet bb D-wave 
states, where r(T(l^L'j) — )• LH) are predicted by us in Eq. (j56ap . and the decay widths of the 
other modes are the same as those in Table IX of Ref.|44] 



Level 


Final state 


Width (keV) 


B (%) 






21.4 


53.1 






11.3 


28.1 




/ 1 AO \ ^ / 


0.58 


1.44 




LH 


6.91 


17.2 




Tvrvr 


0.07 


0.17 






0.0015 


0.0037 




all 


40.3 


100 




/ 1 AO V -L / 


22.0 


77.1 




^ _L y.flSp ) 

/ 1 AO V ^ ^ } 


5.7 


20.0 




LH 


0.75 


2.63 




Ttttt 


0.07 


0.25 




all 


28.5 


100 


T(l3i?3) 


7 + X6(l'P2) 


24.3 


89.6 




LH 


2.75 


10.1 




Ttttt 


0.07 


0.26 




all 


27.1 


100 



B. LH Decay of ^(3770) 



Recently, BES reported 48|-|50| that the branching ratio of the non — DZ) decay of ^(3770) 
is about 15%. While the corresponding data of CLEO 5l|] imply zero. The total width 
r('?/'(3770)) is 23.0 ±2.7 MeV 521"^, and the hadronic and El radiative transitions contribute 
about only 350-400 keV and 1.5-1.8% to the decay width and the branching ratio of non- 



In this subsection, we still cite PDG06 data, to be in consistent with our analysis in Ref. 4C | 
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{ 



} 



FIG. 12: QCD Feynman Diagram for the S-D mixing term. 



DD decay mode, respectively. To clarify this puzzle, the annihilation decay of ^(3770), i.e. 
'?/'(3700) — 7- LH, is considered in our previous paper[40|, where '0(3770) is taken as a D-wave 
dominated state with a small admixture of the 25" state. We found when the annihilation 
decay is included, T{'ip{3770) — )■ non — DD) is 1.15 ~ 1.20MeV, corresponding to branching 
ratio of about 5%. 

In the above sections, the short- distance coefficients and long-distance matrix elements 
of iIj{1^Di) have been given in detail. Now, we show how to get the S — D mixing term. The 
typical Feynman diagram for interference between the color-singlet ^Si and ^Di is shown in 
Fig. [12]. The interference between other Fock states of 5*- wave and D-wave are suppressed 
by as or v"^. For example, the interference between two P-wave octet states is of relative f ^ 
order. And the S-wave singlet interference is of relative order, since there are at least two 
additional gluons in the S-wave Fock state of '^D^^^ In full QCD, the square of the D-wave 
amplitude is logarithm divergent in phase space integration, and that of S-wave amplitude 
is finite, therefore, the combination of them will be finite. Then the short-distance part in 
Eq.( l59|) could be calculated in 4 dimension: 

2lmA{{QQfl^ ^ (QQ)il) = Re[J2 iMiQQth, ^ mM*{{QQ)^l ^ LH)]rf<|.. 

(59) 

Taking into account the corresponding long-distance part, we then obtain the final expression 
for the mixing term in Ref. 40| : 



VI. SUMMARY 

In this paper, in the framework of NRQCD we study the light hadron (LH) decays of 
the spin-triplet (5=1) Z^-wave heavy quarkonia. For completeness, the short-distance coef- 
ficients of all Fock states in the ^Dj{J = 1,2,3) quarkonia including D-wave color- singlet, 
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P-wave color-octet, and S'-wave color-singlet and color-octet are calculated perturbatively 
at order. The infrared divergences in D-wave singlet states are absorbed by the P-wave 
color-octet matrix elements. The operator evolution equations of the four-fermion oper- 
ators are also derived and are used to estimate the numerical values of the long-distance 
matrix elements. We find that for the cc system, the LH decay widths of ipll^Dj) pre- 
dicted by NRQCD is about 2 ~ 3 times larger than the phenomenological potential model 
results, while for the bb system the two theoretical estimations of r{T{l^Dj) — > LH) are in 
coincidence with each other. 
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